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Abstract . After a brief description of 
the modern instrumentation used in Geodesy 
for the detection of the d e forma t i ons of the 
crust of the Earth, follow a list of problems 
and a discussion concerning the needs for 
the survey of the physical quantities of in_ 
terest in geodesy, geology, geophysics and 
engineering such as the strain invariants, 
the optimal network of baselines and the a£ 
curacy. An analytic method is also given for the 
computation of the effect of a source of dilatation 
in a spherical Earth. 

In t roduc t ion . 

From the observation of the stress drop 
in the earth’s crust caused by earthqua- 
kes in several regions of the world we have 
learned that the most frequent stress drops 
are those at low level (Caputo, Console., C1978)); 
but we have also learned that stress drops 
ofless than 0.4 bar are almost never found 
in earthquakes recorded with usual S’eismo- 
graphs. The corresponding strain released 
amounts to several parts of 10“^. This is a 
quantity which is measurable with instrume- 
nts of various type, including portable ge£ 
detic instruments, with an accuracy better 
than 10%. This in turn implies that we can 
try to monitor the accumulation of strain 
which preceedes the earthquakes. The same 
accuracy is sufficient for monitoring the 
deformation of large structures such dams, 
bridges, buildings, tunnels, mineshafts and 
others. To determine this accumulation of 
strain one should observe directly in the 
rock. This in turn can be achieved with se- 
veral techniques some of which are of geo- 
detic interest. Other geodetic techniques 
are of lesser accuracy but are of interest 
because they can be applied in some circum- 
stances when the more accurate techniques 
are not possible because of logistic diffi- 
culties; they will be discussed in a separa- 
te presentation.' 

In the following we shall first describe 
the geodetic techniques used. in the past, 
and discuss for future needs. Some advances, 
will also be presented. 

Geodetic measurements . 

The basic concept in the geodetic meas_u 
rement of horizontal deformation is that 
the distances between selected points are 
measured at different times and then com- 
pared. The instrumentation is based on el- 
ectromagnetic waves which are used in two 
ways . 

One way is to observe the motion of the 
interference fringes of a monochromatic 
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light source caused by the relative motion 
of two bench marks connected with the light 
source and a reflector. This method allows 
to monitor continuously slow variations of 
distances with a resolution of the order 
of 10”^^. The distance between the bench 
marks is limited by the fact that the light 
beem should stay in the vacuum; distances 
of almost 1 km have been used. There should 
not be basic difficulties in making this 
instrument portable or to operate it betw- 
een bench marks connected to the rock such 
as in natural caves in mine shafts, rail- 
way of road tunnels, or in tunnels drilled 
for this purpose. Since the in s t rumen t doe s 
not give any absolute measurement of dist- 
ance but it monitors its variation the 
need to make* it portable is for saving time 
in assembling it and to make its use easier. 
A good example for this type of instruments 
is that of Wyatt and Berger (1978). A me- 
thod of much lesser accuracy, at most one 
part in 10”^, but of much easier and faster 
use, is that of measuring distances between 
several points using laser instruments with 
one br more colors. To increase the accuracy 
of the instruments with one color Slater 
and Huggett (1976) introduced the technique 
of using a radial set of baselines, of al- 
most instantaneous measurements and of com 
paring the ratios of the lengths of the ba 
selines at different times to detect rela- 
tive variations of length. The method redu 
ces the errors due to the variation of the 
atmospheric parameters and is applicable 
with the hypothesis that the longterm chan- 
ges in the index of refraction over the ba- 
selines are proportional to the index of 
refraction at each line; the proportiona- 
lity coefficient must be the same for all 
the lines in the array during the specific 
survey. The method of the ratio is useful 
because it requires to record the atmosph£ 
ric parameters during the survey only at 
one end point. of the baseline. The preci- 
sion of the method is about one or two parts 
in 10“^ . 

This method has been succesfully used 
by measuring the atmospheric parameters at 
both end point of al 1 . base 1 ines for very 
long baselines (of the. order of 100 km); 
the precision of the method is reported of 
2 parts in 10”^((Carter and Vincenty ( 1978)). 
Another method to minimize the errors due 
to the non uniform and constant atmosphere 
has been introduced by Farm (1973, 1975); 
he used . base . 1 ines over flat terrain and 
made the observations only when the wind 
was almost parallel to the baseline. 

These methods do not give the absolute 
measure of the baseline with the same ac- 
curacy obtained measuring the average at- 
mospheric parameters along the baseline by 
means of an aircraft (Savage (1975) and Sa- 
vage and Prescott (1973)) although the va- 
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riations of the length of the baseline is 
obtained perhaps with greater accuracy. 

However the best results in measuring 
the variations of the length of baseline 
and its absolute length are obtained with 
a multiwave length geodimeter (Hugget et al . 
(1977), Slater and Hugget (1976)) with an 
accuracy of 10 ”^ in the variation of the 
length and in the absolute value of the 
length of the baseline. It seems therefore 
that the accuracy of the instrumentation 
would allow to obtain the necessary accu- 
racy on the surface of the earth. 


The computation of the strain invariants. 


Problems which are still partly or com- 
pletely unsolved are those connected .with 
logistics and cost. But the most important 
problem to solve is the meaning of the in- 
formation which., is now being collected on 
the surface. of tHe earth in connection 
with what is happening below it and which 
we really want to know. Strictly speaking 
the problem which I should discuss here is 
two dimensional but in reality the problem 
is intrinsically threedimensional and one 
should never forget it. Therefore. I shall 
try to set it. in its proper space, the pa£ 
tition line between horizo.ntal and verti- 
cal deformations will be set later where 
and how it will seem most proper. 

The stress accumulated in the Earth’s 
crust probably has uniform distribution 
over fairly large scale at least in the 
tectonic processes. This circumstance al- 
lows to observe all. the components of the 
stress and therefore the stress invariants 
even when the . base 1 ine s of the trilatera- 
t ion are . large . . 

If one succeeds in measuring the displa- 
cement vectors 

U 32 ). U 3 (“13 .U 2 3*“33^ >’‘ 2 j ’’'3j^ 

^ £ j are cartesian orthogonal coordinates with 
the axis normal to the surface of the 
Earth;, assuming that the displacement is 
a linear function of the coordinates, one 
may write the following system of 9 equa- 
tions . 

du . . 

u. .-w. = (x, .-X. ) ,i,j,k,=l,2,3 (1) 

J.1 JO dx^ k-i ko 

du^ 

in the 9 unknown which will give the 

components of j the strain and therefo- 
re the invariants. It may be noted that on 
the surface of the Earth one may reasona- 
bly assume 

du« du. du^ du« du^ . du- dn« 

3 _ ^ 1 3 ^ ,2 3 ^ s X . 1 ^ 2 V 

dx^ dx^, * dx^ dx^’ dx^ si -1 nx^ dx^ 

where si is the Poisson ratio. It is obvi- 
ous that if one observes only horizontal 
displacement and one estimates the invari- 


du^ du^ 

ants with the hypothesis that -3 ” -j * 

== 0 , with few exceptional ^ 3 ... 

cases (e.g. see Caputo 1978), the results 
could.be misleading. The cavity effects and 
the topographic and thermal effects (Harri- 
son 1976, Harrison and Herbs t 1977) give 
classic examples of wrong interpretation 
of strain observed in wrong places or bet- 
ter they give examples of wrong extrapola- 
tions of data, as it would be in our case. 

Since we can measure d isp 1 acemen t s on 
the surface of the Earth the question is: 
how may we observe the displacement vector 
in points under the surface of the Eargh? 

An estimate of the vertical component 3 

at depth cpuld be made by measuring dx- 
directly the dilatation d and then subtract 
ing the surface components computed with 
measurements on the, surface of the earth. 

The Sacks Everston (Sacks et al . 1971) di- 
latometer is a good example of an instrum- 
ent which can measure dilatation with an 
accuracy of 10“^. But the additional know- 
ledge of d does not allow to compute the 
maximum shear strain and its direction which 
are of the. greatest importance. The real 
problem is to estimate 


du., du du^ du^ ■ 

diq j^...levellings allow to 

estimate the last two and therefore to retrieve 
the first two, but only on the surface. 

Operating in a mine shaft could be an expe- 
riment to try, however, in this case, the 
displacement at the end of the baseline 
could be determined only by means of an 
open poligon with an accuracy which would 
hardly be acceptable and in any way would 
be by far inferior to that of the surface 
measurements. The use of a road or railway 
tunnel could allow to determine the displa- 
cement vector with a closed poligon, but in 
most cases the tunnel would not be located 
in a useful position. The problem is open 
and it is one of the most important. 

A problem which I think we should also 
consider is that of the most suitable di- 
stribution of the baselines or of the bench 
marks for determining the distribution of 
strains on the surface of the Earth in the 
case of the hypothesis du,/dx^ = du /dx^= 0 . 
In this case the system (1) is reduced "^to 
2 linear systems of 2 equations each namely 


"11 

'"lo 

dUi 

dUi 

dx7^*2l"’'20^ 



du 

dx7^’'22“*20^ 

“12- 

""lO 


“21" 

''20 = 

d“2 

dxi 

^^12 dX 2^’'21 * 20 ) 



du. 

du2 

“22“ 

^20" 


^*12~*10^* d7J^*22"*20^ 
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The matrix of the coefficients in both sy- 
stems is 


11 

^10 

^21 ' ^20 

r 

(3) 

12 

^10 

^22 “ ""20 



which should be different from zero.; The r£ 
fore the three points pQ,P^,P 2 /"should noX 
beinastraightline. 


Accuracy of the d e t ermina t ion o f strain , 


Given the m.s.e.s (x^)=s(x 2 )in the deter_ 
mination of the coordinates of the points 
of the network in all the surveys, which 
in turn give the displacement vectors, th'e 
m.s.e. of the i nvar i an t s ^ ar e in the first 
approxima tion 


s(d) 


s 

sin z 





2s 


J 


s(a) = s 
s ( S ) = 2 V7^ s 


(4) 


Where S is the max imum - shear strain, a is.- 
its direction, z is the angle P^ P^ and 

d, and d are the length of P, and P^ P_. 
If the strain field is uniform then the best 
s(a) and s(S) are obtained for z = tt/ 2 and for 
dj and d^ very large. This implies that when 
one computes the strain invariants one loses 
at least a factor 3 in the accuracy. 



Fig,, la - Dilatation as function of time in 
the period 1970-1977 in the Phl£ 
grean Fields near Naples. Each 
diagram gives the dilatation in 
the point connected by the dash 
ed arrow. In the top left diagram 
are indicated the units of dil£ 
ration and the time interval in 
which it has been accumu^lated . 



in the period 1970-1977 in the 
Phlegrean Fields near Naples, 
Each diagram gives the maximum 
shear in the point connected by 
the dashed arrow. In the top 
left diagram are indicated the 
units of dilatation and the time 
interval in which it has been a£ 
cumu la'te'd',' 



Fig. Ic - Network of baselines for the me£ 
surement of displacements across 
the Messina Straits. 


When the strain field is not uniform 
formulae (4) set a limit to the accuracy of the 
determination of the invariants. It is the 
indetermination principle of the strain fields 
observations, which is expressed by 
^ e > s , 

where ^ is the strain observed in the path - 
An almost complete example of a determination 
of a strain field is that of the Phlegrean Fields 
near Naples (Caputo 1978) shown in Fig. la, lb. 
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Fig. Id - Displacement reported in the two 
points located North in the net- 
work of Fig . 1 . 
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Network of baselines for the mea 
suremnt of strain in the Ancona 


The accuracy which one may achieve with 
standard portable instrumentation and pro- 
cedures, in a triangular network with 5 km 
baselines is s 1 cin and one ob'tains s(d) = 
=6*10’6; if one uses the mekometer then 
s ^ 2 mm and the baselines are 2 km there- 
fore s(d) = 3‘10”^; by using a multiwave- 
length gepd ime t er with baselines of 5 km 
s(d)’would be 3 ‘10“^. 

It is desiderable to develop numerical 
methods and network designs in order to im 
prove the accuracy in the determination of 
the invariants. 


Fig. If - Network of baselines for the me 
suremnts of displacements betwei 
the Appennines and the Alps. 


With reference to the network design 
let us consider the limiting case of a ne 
work of infinite size in two dimension. I 
the distribution of the strain is unknown 
obviously one should consider the most sy 
metric distribution of bench marks. Accor 
ding to the principle of the Pflastersatz 
the most simple distributions of points 
which have properties of symmetry are the 
three indicated in the Fig. 2 

When one considers the baselines neede 
to determine . the invariants of strain one 
sees that the networks (c) of Fig. 2 is no 
convenient because of the length of the b 
selines which should also be VT” 1 (with 1 
the minimum distance between points). 

Since the two networks of points (a) a 
(b) have no. condition equation the ratio 
of the numbers of measured baselines to 
the number of points where the strain is 
measured is the same. Also, to the first 
order s(d) = for the network (b),^3 

s(d)= M for (a)^0; that is they are almost 
the same. However the distribution of the 
points where the invariatns are measured 
is better in (a), because in (b) some poi 
nts have the invariants which are direct 
interpolation of others. 

We may also note that it may be needed 
to reinforce the strength of the networks 
which we have considered here and for whi 
ch, it is immediately seen that the condi 
tion equations are non existing. If the ' 
baselines are short enough one may consi- 
der to measure some of the angles or some 



(C) 

Fig. 2 - Possible networks of baselines 

formed with regular poligons. 

more baselines between the existing points. 
For instance in the network which we used 
across the Messina Straits (Caputo et al. 
(1974)), when we measured only the baseli- 
nes we had only one condition equation; 
we had more condition equations when (at 
times) we measured also the angles; the 
decision to adopt this procedure was due 
to logistic problems. In the case of the 
network of the Ancona area (Baldi et al . 
(1976)) we decided instead to add more poi- 
nts to the network and we had 9 condition 
equation. The method of adding more base- 
lines is often preferred to the measure- 
ment of angles because one obtains the 
strains in a larger number of points and 
therefore obtain a higer resolution in the 
distribution of strain; also with the mo- 
dern instrumentation, in networks with ba- 
selines of 10 km or longer, the measurem- 
ent of angles is almost one order of magn^i 
tude less accurate. Before closing I must 
call the attention to the set of problems 
which I raised and which we should eventu 
ally discuss with the many more which pr£ 
bably will be raised after my presentation.. 
These problems are 

a) the need to have a threedimensional 
description of the strain fields 

b) the need to have available a multi- 
wave length geodimeter 

c) the need to monitor strain on long 
baselines with an interferometer ea- 
sy. to install and to operate 

d) the need. to measure on bench marks 
which are rigidly connected to the 
medium to observe 

e) to establish networks which are opti- 
mal for the final objective 

f) to select a test ground where to make 
an expe r imen t to satisfy. all these r^ 
quirement and to do the experiment. 
(This is a proposal for those who are 
interested). 


g) to develop mathematical models for the 
displacement fields caused. by internal 
perturbation of the Earth elastic 
field. 

h) the use of theodolites, in conjunction with 
electromagnetic distance equipment with 
accuracy s, is of use only if the length of the 
baseline 1 is 1 < 8 * 10 s in which case its 
error would be 2 s. 

i) preference to interferometers when investi- 
gating slow rate accumulation of strain, 

j) need to adjust the observed displacement 
fields with the Navier equations. 

The displace rent field caused by a sour ce 
of dilatation in a spherical Earth . 

It has been established that pryor to 
some earthquakes there. is a deformation 
of the. Earth’s crust measurable on its 
surface. According to th e ^ d i 1 a t ancy model 
in some cases this deformations caused by 
a dilatation in a threedimensional. port- 
ion of the Earth’s crust. It is therefore 
needed to have a threedimensional geome- 
tric elastic model of this phenomenon for 
the - interpretation of the geodetic data. 

A solution, for a flat Earth model and 
a point source of dilatation. has been g^ 
ven by Sacks et al (1975); we shall ex- 
tend here . this solution to a spherical 
Earth model S of radius r and a point 
source of dilatation at a distance h 
from the center of the Earth. 

The method. consists in finding first 
the field u generated by the dilatation 
source in P in an infinite space assum- 
ing the origin. of the polar coordinate sy- 
stem in P^; then. we shall move . the origin 
of the coordinates in the center 0 of S 
and finally add. to u a vector v which ren 
ders the component of the stress applied 
to S nil on S as it naturally is. In a 
mo r e ° s o f i s t i c a? ed solution one may. render 
t such to balance the weight caused by 
tfie topography. 

We shall assume that the body forces 
caused by the. source and its perturbation 
are negligible; in the first approximat- 
ion this is acceptable when the source 

dimensions and the dilatation are limited 
to realistic cases. We assume also that 
the gravity field of the Earth has negli- 
gible effect with. respect to . the other 
fields playing the major role in this phe 
nomenon. Solutions taking into account " 
all those effetcs would obviously be more 
accurate; the solution presented here is 
considered only a first approximation of- 
the spherical Earth. 

The displacement caused in the space 
by a generic source if dilatation in P 
with the condition of co n v ergence at infi 
nity is “ 

u^= dV r = 0 - (5 ) 
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where dV is the change in volume of any 
volume which includes P. . It is immediate 
ly verified that the dilatation. of this 
displacement field is nil at all points 
different from P^. 

. The formule for the mentioned transfo£ 
matioh of coordinates , with reference to 
thefigureS are 

t [ 2 .2 ^ h sin z 

r = r +h -2rh cos 2 r = — : (6) 

L J sxn a 


1 and m are the elastic parameters. 

The components of the stress of the s£ 
lution V. are 


n |j^(n+3)l+2(n+i)mj ^-ln(n+l.)B^ | r"+ 


+ 1 (n) l+2(n-l)m 


1 [ 


-n(n+l)l J.r 
J 3,n 3,n^ 


n-2 


Pn 


and substituting in (I) we obtain 

[ 22. 1-i 

r +h -Zirh cos z cos a 

(7) 

Uz=dV ^r^+h^-2 r h cos ^ J ^ ^ 

where sin a and cos a are obtained from 

(6) : 


r 2 2 * T -i 

sin a = h sin 2 r +h -2 r h cos z 




• 2 

h sin 2 


( 8 ) 


2 2 

+h -2 rh cos z 


dV(r - h cos z) 
j^r^+h^-2r h cos 

dV h sin z 
^ |^r^+h-^-2r h c 


( 9 ) 


■] 


3/2 


* (31 A^.O* 2 

m Bj^ 

00 


t = m Z] n 

(n B 

rz , 1 

A 1 n-2l 

dPn 

""3,nJ J 

dz 


CD 

+ 



Fig. 3 “ Coordinate system used in formu 

lae (5) through (f2) • 


We shall compute now the solution v of 
the Navier equations which is to be added 
to u in order to obtain the. final solut- 
ion u + V such that the stress components 
t applied to are nil on S . 

According to the formulae of 8aputo (1961) 
the general solutions of the Navier equa^ 
ions in the axisymmetrical case and homo- 
geneous Earth of this case are 




\ 

+A 


n-1. 2n+i _ . . 

r +A. ^ r ) — P (cos z) + 

3 ,n z n 


(10) 


iO 


00 

■E„(b 

I " l>n 


r +A 


3 ,n 


n-i. 2n-H 
^ ^2n(n+i) 


B 


3,n 



dPn 

dz 


B 


i,0 


3g-^5 

2 



If we require that the components of the 
stress t of the solution u + v applied 

to S are nil on S the conditions are 
o o 


In 

2 


I ^(n+3) i+2m(n+i) J A^ ^-ln(n+i) B^^ 


n-2 


(n-i)(l+2m)A - 1 n(n+i) B» I ro | Pn 

3,n 3,nJ 


+ 3 1 A 


+ 

’r*=ro 




ro sin z dz 


y' [(A, +nB )r^ +(A^ +(n-2)B^ )r ^ ^1> 

^n I l,n i,n o 3,n 3 ,n o 
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dPn [duz _ uz . 1 ^1 _ „ 

"^~dT * ["dT T 7 dz J ■ ° 

r~r 

o 

where A. and B are the unknown. 

. 1 , n . ij,n 

An expansion of u and u , appearing 
in (1^, in series of Legendre polinoraials 
gives finally the values of A. and B . 
to sobstitute in (10) to obtain* ?he finVl 
so lu t ion . 

The extension of this result to a pre- 
stressed, Earth is almost obvious. What 
we really need is to prepare tables and 
figures with the computations of the 
effects of a dilatation in according 
to the formulae and to observe the phen£ 
mena on the real Earth. 

The extension of the results of this section 
to a layered Earth is also trivial (Caputo 1961). 



Fig. 4 - The four points P, should not belong to a 
single plane. 
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